Parity detection and entanglement with a Mach-Zehnder interferometer 
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A parity meter projects the state of two qubits onto two subspaces with different parities, the 
states in each parity class being indistinguishable. It has application in quantum information for 
its entanglement properties. In our work we consider the electronic Mach-Zehnder interferometer 
(MZI) coupled capacitively to two double quantum dots (DQDs), one on each arm of the MZI. 
These charge qubits couple linearly to the charge in the arms of the MZI. A key advantage of an 
MZI is that the qubits are well separated in distance so that mutual interaction between them is 
avoided. Assuming equal coupling between both DQDs and the arms and the same bias for each 
DQD, this setup usually detects three different currents, one for the odd states and two for each 
even state. Controlling the magnetic flux of the MZI, we can operate the MZI as a parity meter: 
only two currents are measured at the output, one for each parity class. In this configuration, the 
MZI acts as an ideal detector, its Heisenberg efficiency being maximal. Initially unentangled DQDs 
become entangled through the parity measurement process with probability one and for a class of 
initial states our parity meter deterministically generates Bell states. 

PACS numbers: 73.23.-b, 03.65.Ta, 03.67.Bg, 03.65.Yz 



I. INTRODUCTION 

A parity meter couples two qubits and allows to mea- 
sure whether the qubits are in the even subspace spanned 
by the even states {|tt)J4-i)} or in the odd subspace 
spanned by the odd states {Iti), lit)}- A particular ba- 
sis for these two subspaces are the four Bell states, which 
are the maximal entangled states for two qubits. As these 
Bell states are also eigenstates of the parity operator for 
two qubits, a parity measurement drives the qubits into a 
Bell state and constitutes a source of entanglement. The 
parity detection is thus of special importance in the possi- 
ble implementation of quantum computation schemes ii^— 

Realizations of parity meters for charge and spin qubits 
have been discussed considering a quantum point contact 
(QPC) as detector )^i^ The purpose of our work is to inves- 
tigate parity detection of charge qubits with the help of a 
Mach-Zehnder interferometer (MZI) as shown in Fig. [T] 
The MZI is the simplest possible type of interferometer. 
It has recently been implemented in a two-dimensional 
electron gas in the quantum Hall regime with electron 
motion along edge states<^"— A discussion which com- 
pares the MZI with a second order two-particle interfer- 
ometer is presented in Ref. Il2l . The non-linear behavior 
of the MZI has received much experimental^"— and the- 
oretical attentionJ^"— In the following we assume that 
the applied voltages are small and the MZI is operated 
within its liirear response regimeJi The MZI coupled to 
qubits in its arms has the advantage that the two qubits 
are (on the mesoscopic scale) far away from each other. 
Therefore they are not subject to mutual Coulomb in- 
teractioir. An additional advantage is that the MZI is 
tunable both through a variation of the applied magnetic 
field and through gate voltages. These gate voltages al- 
low to determine the length of the interferometer arms 
and to control the transparency of the QPCs which act 
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FIG. 1: Mach-Zehnder interferometer (MZI) coupled capac- 
itively to two qubits. xi(X2) is the phase accumulated by 
the electrons passing through the upper (lower) arm of the 
MZI over a length L. These phases depend on the state of the 
qubits. C\ , Ci and Ci are the capacitances characterizing the 
double quantum dots (DQDs). 



as beam splitters. Admittedly the fabrication of such a 
sample is not simple since one of the qubits would have 
to be built in the interior of the Corbino disk. As a con- 
sequence its gates can only be biased with wires which 
will bridge the conducting region of the interferometer. 
Our analysis of the parity detection is based on a for- 



mulation which takes into account that a quantum mea- 
surement lasts a finite time. This assumes that the de- 
tector is only weakly coupled to the system. Such weakly 
coupled detectors are typical for on chip detectors in 
mesoscopic physics and have led to a strong theoreti- 
cal development ji^"— Previous works have shown that 
linearly or quadratically weakly coupled detectors such 
as QPCs distinguish in general three currents. Their 
measurement process allow purification of a mixed state 
and entanglement of the qubitsi^^— In subsequent work 
Trauzettel et al^ have shown how to use a QPC as a 
parity detector. The investigation of entanglement gen- 
eration finds that the detector still distinguishes three 
states^, two states which are detectable in repeated mea- 
surements belonging to the two parity classes and a 
time-dependent state which oscillates between the parity 
classes. These results follow from a qubit Hamiltonian 
formulated in the iTj^-basis. In this work we treat qubits 
without tunneling and show that the detector can be set 
to measure only two currents. 

The paper is organized as follows: In Section II, we de- 
scribe our setup (see Fig. [T]). In Section III, we give the 
conditions for the MZI to operate as parity meter and 
in Section IV, we investigate its Heisenberg efficiency. 
The efficiency is determined by the ratio of the mea- 
surement time of the detector and the dephasing time 
of the qubits i^i2^i2^ We find that a properly tuned MZI 
parity detector has maximal efficiency, i.e. it is a quan- 
tum limited detector. In Section V, we show that the 
measurement performed by the MZI as parity meter en- 
tangles the DQDs with probability one when each DQD 
is initially in a superposition state. We also investigate 
in this last section another interesting property of this 
parity detection: for this particular class of initial states, 
if the odd current is measured, we know into which odd 
Bell state the qubits are driven. This makes the MZI an 
efficient parity detector and, for initial superposed states, 
a deterministic source of entanglement. 



II. MODEL 

The electronic MZI is built in a quantum Hall bar by 
using two QPCs as beam splitters. The qubits are spa- 
tially separated so that they do not interact directly with 
each other. The interaction takes place between each 
qubit and the charges of the corresponding arm and is 
described with the help of geometrical capacitance coef- 
ficients, see Fig. [T] Ci determines the coupling of the 
charge on the DQD to the charge on the arm, Ci deter- 
mines the coupling between the dots, and C2 is the capac- 
itance of the DQD to the ground. For simplicity, here, 
the capacitances are assumed to be the same for each 
qubit. The regime of weak coupling between the detector 
and the DQDs is defined by Ci <C C2,Ci. Additionally, 
to control the system, we can tune the Aharonov-Bohm 
flux, the length of the arms, and the transmission and 
reflection probabilities of the QPCs. The Hamiltonian of 



this model reads 



H — i/qb + ^int + ffdot, 



where 



H, 



qb 



<^1 -1 ,'^2.2, ^1 ;^1 , ^2 .2 
.l + 'T^ 



^int = Ul{^-^)+U2C-^). 



(1) 

(2) 
(3) 



Hdct is the Hamiltonian of the MZI. _ffqb is the Hamil- 
tonian of the qubits. We assume no tunneling between 
the dots, Ai = A2 = 0, and the same bias for the two 
DQDs, £1 = £2 = £, so that iJqb = (e/2)(a-] + &l). The 
Pauli matrices for the upper and lower qubits are a\ and 
(7^. In the interaction Hamiltonian Hmt the expressions 
(1 — CT^)/2 and (1 + (T^)/2 are the charge operators of 
the qubits which couple to the charge on the arms of the 
interferometer. The Mach-Zelmder detector considered 
here is a linear detector in the sense that qubits couple 
linearly to the charge in the arms of the MZI through the 
dz operator. As the qubit Hamiltonian is also in the a^- 
basis, there is no relaxation process for the DQDs due to 
the coupling with the detector. As shown in more detail 
in the next section investigating the MZI as parity meter, 
the parity operator P = a\®a1 commutes with both i?qb 
and iJint- It follows that, in this conflguration, the Mach- 
Zehnder detector is a quantum non demolition detector 
(QND-detector) of the parity classesi^"— This is to say 
that repeated measurements will find the qubits in the 
same parity class, the MZI detects only two stationary 
currents. This is an important distinction between previ- 
ous works on entanglement generation which use a qubit 
Hamiltonian in the (Ta;-basi9^i2ir— and the investigation 
presented here. To precise ifint, U1.2 are the potential op- 
erators on the upper and lower arms respectively. They 
depend on the bare charge and on the screening charge 
on the arms due to the presence of the DQDsi^i^ This 
interaction Hamiltonian implies that the electrons tran- 
siting in the arms of the MZI are sensitive to the charge 
in the DQDs only if this charge is in the dot close to the 
arm. As a consequence, electrons on arm i = 1,2 will 
acquire during their transfer a phase Xi due to the mag- 
netic flux, due to the length of the arm i and due to the 
state of the DQDs: 



Xi = fc(d - L) + 6'„p + Axi 
X2 = fc(d - i) + 6*^ + Ax2 



2 

1 + ^2 



(4) 
(5) 



Axi is the additional phase depending on the state of 
the DQDs and 6up — Od = 27r$/$o is the magnetic flux 
between the arms. For simplicity we consider both arms 
having the same length d and the electrons are sensitive 
to the DQDs over a length L, k is the Fermi wavevec- 
tor, see Fig. [TJ Equations ^ and ([5]) should be viewed 
only as abbreviations, since in reality the phases and the 



transmission probability T31 to be given in the next sec- 
tion are c-numbers depending on the state of the qubits. 
Thus instead of writing operators in Eqs. H} and ^ 
we should index these quantities with the corresponding 
qubit states. Here we avoid such a heavy notation and 
keep operators. 

In the next sections, we will derive the conditions un- 
der which the MZI coupled to DQDs allows their parity 
detection and operates as ideal detector in this regime. 
We will then investigate the creation of entanglement 
through the parity measurement. 



be in general distinguishable: one for the odd states lo 
and two for the even states I-^-^ and I^ as shown by the 
three black points on Fig. [2]a. Although this measure- 
ment is not a true parity measurement, it already allows 
entanglement between the qubits as investigated with a 
QPC coupled to two DQDsj^i^i^^ In our work we want a 
true parity meter with advantages for the entanglement 
probabilities as discussed in the next section. Equations 
© and (|12p and Fig. [5]b show that measuring only one 
current If, for the even states is possible if the magnetic 
flux is set to 



III. PARITY METER 

A parity meter is a quantum detector which can only 
distinguish states belonging to different parity classes, 
the odd and even classes. The observable we measure is 
the averaged current in lead 3 (see Fig. [1]). Assuming 
that the electrons are injected in lead 1, the current is 
given byi^ 



27r$/$o = modulo tt. 



(14) 



2eV 



(6) 



T31 is the transmission probability from lead 1 to lead 3, 



T31 = {TlRr + TrRl + 2^ TrRlTlRr cos $) . (7) 

It can be calculated with the help of the scattering matrix 
S of this model: T31 = 15*311^, where 5*31 is one element 
of the S'-matrix. It corresponds to the amplitude proba- 
bility of going from lead 1 to lead 3, 



^3 



(8) 



We introduce the phase $ = Xi — X2 where xi and X2 are 
given in Eqs. ^ and ([S]). Depending on the four possible 

states of the qubits {|tt>, Itl), lit), III)}, ^ can take four 
different values: 



$tt 


= 27r$/$o-Ax2, 




(9) 


*n 


= 2^$/$o, 




(10) 


*it 


= 27r$/$o + Axi - 


-Ax2, 


(11) 


*u 


= 27r$/$o + Axi. 




(12) 



To operate the MZI as parity meter, only two currents 
have to be distinguishable: Iq corresponding to the odd 
states { Iti) , lit) } and le corresponding to the even states 
{ltt),IU)}- Considering Eqs. (jlOlllip . only one current 
lo for the odd states will be measurable if the incremental 
phases Axi for electrons in upper and lower arms are the 
same: 



Axi = Ax2 = Ax. 



(13) 



Assuming Eq. ()13p and measuring the output current 
as a function of the magnetic flux, three currents will 



providing a maximum or minimal transmission. To be 
specific, in the reminder of the work, we will choose 
27r$/<I'o = 27r to be at the maximum of transmission. 




27r$/$c 



b) 




27r$/$o 



FIG. 2: Current measured at lead 3 as a function of the mag- 
netic flux between the arms of the MZI. a) At an arbitrary 
flux three currents /o,/tt ^'I'l ^ii ^''® distinguishable even 
if the incremental phases are identical: Axi = Ax2 = Ax. 
b) To measure only two currents corresponding to the parity 
classes, lo and le, we set the flux to 27r'l>/<I>o = modulo vr. 
A maximum of the current is reached at 27r$/<l>o = 27r. 

Equations ([T^ and P^ are the two conditions re- 
quired for the MZI to act as parity meter. Then the 
transmission probability T31 can be expanded around 
27r$/$Q = 27r up to second order in Ax. It depends 
then on the states of the qubits only through the parity 



operator P — a\ ® al: 
^31 = TlRr + TjiB-L 



+2^TrRlTlRr{i - -^il + al(g)al)y (15) 

This measurement operator P allows parity measure- 
ment and the creation of entanglement. Indeed P has 
two eigenvalues ±1 of degeneracy 2, corresponding re- 
spectively to the even parity class and to the odd parity 
class. Its four eigenstates are the Bell states, 

i^+) = ±(\n)+\ii)), iVe")=^(in>-iu»,(i6) 
i^o+> = -^(\n)+\m, iV'o->=^(in>-i;t».(i7) 

The lower indices indicate the parity class and the upper 
indices (±) indicate whether the sum or difference of the 
computational states is taken. These Bell states are of 
particular interest as they share the properties to be the 
maximal entangled states for two qubits and to have a 
definite parity. A key property of the qubit Hamiltonian 
-ffqb considered in this work is that its dynamics conserves 
the parity of the Bell states: 



H^kHt)=^H'f) and iJqbl^^ 



^\^t 



(18) 



It is the parity measurement in the (Jz-basis combined 
with the qubit Hamiltonian Hqh and the coupling in this 
(Tz-basis (no relaxation process) that allows the MZI to 
act as parity meter for the DQDs during the entire mea- 
surement process. In this sense, this weak measurement 
procedure is the analog of a QND measurement of parity 
classes. 

In the next section we will investigate the efficiency of 
the MZI as parity meter to show that it acts as an ideal 
detector. The efficiency rj is defined as the ratio of the 
dephasing time over the measurement time, < 77 < 1. 
If 77 = 1, the quantum detector is an ideal detector: it 
extracts information on the quantum system as fast as it 
dephases it. In other words the measurement process is 
the only dephasing source. The ideality of the detector 
will constitute a key property for the success of entangle- 
ment generation (Section V). 



IV. EFFICIENCY 

The quantum efficiency 77 is defined as 77 = r,„/rdcph , 
where Fdcph is the dephasing (or decoherence) rate of the 
measured quantum system and Fm is the measurement 
rate of the detector. The most efficient quantum detector 
one can implement is a quantum-limited detector, r; = 1; 
it acquires the information about the measured system 
as fast as this quantum system loses its quantum coher- 
ence. Fdeph is the inverse of the coherence time T2 of 
the measured quantum system and Fm is defined as the 



inverse of the time needed by the detector to distinguish 
the signal from the output noise j^ 



F„, = 



(All 
45 



// 



(19) 



Sji is the low-frequency noise spectrun*^ and A/ = 
lo — le is the difference between the even and odd av- 
eraged currents measured in lead 3 (see Fig. [2]b). As- 
suming linear response. A/ <^ Io,Ie, the shot noise is 
independent of the qubits' states: 



"Je + ^o 



2e^V 



->ii 



r31,e,o(l - r3i,e,o). (20) 



Tai e.o are the transmission probabilities for the even and 
odd configurations respectively. They are given by Eq. 
([7]) in which the parity operator P takes the values ±1 
depending on the parity class. The current difference up 
to the second order in the small parameter A^ is 



A/ 



2e'V 



VRrTrRlTl{Ax)^. 



(21) 



The maximal value of Fm is reached for Rl = Rr =■ 1/2: 

r»<«-«-V2)4fA.= 4f (I)', (22, 

Here the additional phase A^ is expressed as a function 
of the capacitances of the system. Ax = C^/C\ where 
C~^ = C^^ +Ci^ +€2"^ + C~^ is the electrochemical ca- 
pacitance of each DQD. Cq is the quantum capacitance, 
it depends on the density of states in the arms of the 
MZIj^i^ This will allow us to compare this measure- 
ment rate to the dephasing rate of the qubits. 
In general, for two qubits, there are six dephas- 
ing rates describing the coherence between the states 
(Itt), Iti), lit), lii)}- One way to calculate them is to ex- 
tract their expressions from the off-diagonal terms of the 
density matrixi^"— i^ Another way consists in studying 
the phase fluctuations of the DQDs arising from the fluc- 
tuating potentials in the armsj^i^^^i^^ Both approaches 
lead to the same rates and the details of the two cal- 
culations are given in the appendix. Sections A and B. 
Taking into account the parity assumptions, Eqs. (J13p 
and (HU, the six rates reduce to three rates between the 
parity classes: 



Feo = 


1 eV 




r 

^ ee " 


eV 
= 2RlTl — 
h 




r 


= 0. 





(23) 

(24) 

(25) 

Tij is the dephasing rate between the parity classes i and 
j, where i and j are either e for the even class or o for the 
odd class. These results reflect the perfect anticorrelation 
between the events on the upper and lower arms of the 



MZI due to current conservation in the system;^ see also 
Appendix B. 

As we are investigating the MZI as parity meter, the only 
dephasing rate of interest for the efficiency of the detector 
is the dephasing rate between the even and odd classes 
Teo- For simplicity we do not give the general expression 
of the efficiency but we only specify it for the particular 
case of Rl = Rr — 1/2 for which the measurement rate 
is maximal. For half-transmitting QPCs, the efficiency 
becomes maximal: 

''"^^'^^ r,o(i?i = i/2) '■ ^'^^ 

The MZI coupled to two DQDs can thus operate at 
the quantum limit if both QPCs are half-transmitting, 
Rl = Rr = 1/2. As shown in the next section, the ide- 
ality of the MZI will reduce the possible sources of deco- 
herence during the measurement process, enhancing thus 
the probabilities to get entangled states. These probabil- 
ities will be extracted from the density operator, calcu- 
lated by solving a stochastic differential equation for the 
evolution of the qubits. 



V. CREATION OF ENTANGLEMENT 

The entanglement of the DQDs through the mea- 
surement process is investigated by solving the differen- 
tial stochastic equation governing the evolution of the 
DQDs^^^ 



Pij 



-[^qb,P]y 



V 45// 

4 



+ ltj]Pi] 



+CW(/. + /.-2^p,,/,)|^. 



k=l 



(27) 



The current in lead 3 is given by I[t) = ^j. Pkklk + S,{t), 
where pkk is the density matrix's element associated with 
the state \ipk)- The indices fc = 1, 2, 3, 4 correspond to the 
four Bell states {\ipt)A''Pe)A''Po)A''Po)} defined in Eqs. 
(HH]) and (fT7)) . Working with the Bell states presents 
a strong advantage when investigating parity measure- 
ments and entanglement as shown in the previous section. 
li corresponds to the average current the state [ipi) would 
produce. As the MZI acts as parity meter, Ii = I2 = le 
and /a = /4 = /o, where I^ and Iq are the even and 
odd currents, ^(t) is the random white shot noise of the 
output characterized by (C(i)C(O)) ~ 2SiiS{t), where Sn 
is the low-frequency noise spectrum of the MZI (see Eq. 

The differential equation Eq. ^27} can be decomposed 
into three terms. The first term describes the standard 
Schrodinger evolution due to the qubit Hamiltonian ifqb. 
The possible sources of decoherence for the qubits are 
considered in the second term of this equation. They con- 
sist of the term (J^ — /j)^/(4S'//) (compared Eq. ([TO|) 1 and 



a term -fij which is finite only if the detector is not ideal. 
As shown in the previous section, the electronic MZI as 
parity meter is an ideal detector so 7^ = 0. The only 
source of decoherence is thus the measurement process. 
The last term describes the backaction of the measure- 
ment result on the qubits evolution through the noise ^(i) 
of the current. In this work, we are not interested in the 
entanglement genesis^ so we look at the time-averaged 
corresponding equation: 



'm 1 ^^> 



I,)' 



45* 



// 



-Pij- 



(28) 



This equation can be solved analytically. In our setup, 
considering a qubit Hamiltonian without tunneling, a 
classical initial state (a fully mixed state) will remain 
classical. Therefore we investigate an initial state in 
which the individual qubits are in superposition states. 



\^,,)^{a\t)+m)®hn) + s\i)), 



(29) 



where q;,/3,7, 5 are complex numbers such that |ap -t- 
|/3|2 = |-)/|2 ^ |^|2 _ I gy definition, the diagonal ele- 
ments of the density operator pa (t) determine the prob- 
abilities to be in the state \ipi) at time t. Thus the prob- 
ability Pi{t) to be in the Bell state \ipi) will be given by 
Pa (t) , the probability to be in the even subspace is given 
by Pe{t) = Pii{t) + P22{t) and the probability to be in the 
odd subspace is given by Po{t) ~ Pssit) + P44(t). These 
probabilities depend in general on the initial state of the 
two qubits given by Eq. (|^^ . From the analytical so- 
lutions of Eq. (|28|) , we find the probabilities to measure 
outcomes in the different parity subspaces: 






(30) 
(31) 



As expected from the qubit Hamiltonian that does not 
mix the parity classes, Eqs. (|30|) and ([3T|l do not de- 
pend on time. The electronic MZI in this setup works as 
perfect parity meter during the entire measurement pro- 
cess. Each time an electron transits through the arms, 
the DQDs are weakly projected onto a subspace with a 
definite parity. In our QND-limit (see Section HI) , we can 
approximate the weak measurement with many projec- 
tive measurements as the backaction does not take place. 
The probabilities for the qubits to be in one of the Bell 
states are given by: 



5+m - l"^ + ^'^l^„„.,2 



p:it) = 



et 
COS" 1 7- 



\ai-m' . 



et 



\ai-m' 



sin^l-p) 
COS' I — 133) 



p: = 

p- = 



\aS+l3-f\^ 
2 



(34) 
(35) 



In the odd subspace the probabilities to entangle the 
qubits depend only on the initial state via the complex 



coefficients a, /?, 7 and 5. This result is again expected 
from the qubits' dynamics which induces no evolution in 
the odd subspace, the odd Bell states being eigenstates 
of iJqb (see Eq. ([TS])). As a consequence, by choosing 
specific complex coefficients for the initial state, one can 
know with probability one in which odd Bell state the 
qubits are if the odd current is measured. The evolution 
is more fussy in the even subspace. Although the qubit 
Hamiltonian conserves parity, it induces an evolution be- 
tween the two even Bell states due to the bias e: 









(36) 

(37) 



If the measurement outcome is the even current, the 
qubits will remain in the even parity subspace during 
the entire measurement process but they will undergo 
oscillations between the even Bell states with a frequency 
e/fi, whatever the complex coefficients a, /3, 7, 5 are. It 
is not possible to determine in which even Bell state 
the DQDs are entangled, they will be in a superposition 
state of 1?/;+) and IV^^)- 

As an illustration of these properties, we consider two 
product states as initial states defined by a,/3,7, (5 = 
±l/-\/2 and we will compare their probabilities. These 
states are not only interesting theoretically^ but they 
also present the strong advantage to be implementable 
with the experimental state-of-the-art of solid-state 
quantum dots.— The first state corresponds to a = (i ~ 
7 = 5 = l/\/2 whereas the second one corresponds to 
a = ^ = 7 = -(5 = l/x/2: 

l^i) - ^(lt> + i;»®(lt) + i;», (38) 

i*2> = ^(it> + i;»®(it>-i;» (39) 

It is straightforward to see that Pe = Po = 1/2 in both 
cases. As expected the even probabilities exhibit an os- 
cillation behaviour between the two even Bell states. The 
probabilities for j^I'i) are: 



pnt) = ^cos^(f). 



,et. 



p-it) = ^sin^J), 



(40) 
(41) 



and the probabilities for |'I'2) are almost the same, they 
just differ by a 7r/2-phase. In contrast the probabilities 
in the odd subspace are time independent and allow to 
distinguish between the two odd Bell states. Starting 
with 1^1), one finds: 



P+ = l/2, p-=0, 
whereas starting with 1^2), one gets: 
P+=0, P- = l/2. 



(42) 



(43) 



This example clearly shows that the electronic MZI act- 
ing as parity meter entangles the qubits with probability 
one (Pg + Po ~ 1) and allows to know exactly in which 
odd Bell state the qubits are if the current lo is measured. 
In a more general way, if one writes the initial coherent 
superposed state for both qubits as 

l*i„> = ^(lt)+e"'M;»®(lt>+e*''i;», (44) 

where 771 and r/2 are phases specific to the upper and lower 
qubit respectively, then the class of states that allows 
deterministic generation of odd Bell states is defined by: 



-r?i =0 



m-m 



3dulo 



(45) 



This property makes it a very interesting quantum de- 
tector and entangler for further experiments. 



VI. CONCLUSION 

With this work, we present an investigation of the MZI 
as quantum detector of two DQDs characterized by a 
bias and no tunneling between the states. Thanks to the 
control of the flux and of the transparency of the QPCs 
that can be tuned to specific values, the MZI operates as 
ideal detector, efficient parity meter and quantum entan- 
gler for the DQDs. Considering the qubits in an initial 
superposed state, the measurement process drives them 
into the Bell states. These states are of great interest 
for quantum information as they are the maximal entan- 
gled states for two qubits. This set of results encourages 
the implementation of such a setup in mesoscopic physics 
but also in other promising fields for quantum informa- 
tion such as circuit quantum electrodynamics. 
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VII. APPENDIX: DEPHASING RATES 

For calculating the dephasing rates given in Eqs. 
I25p . two approaches are possible: the entanglement ap- 
proach (Section A) or the charge fiuctuations approach 
(Section B). The entanglement approach^Si^i^ aims at 
writing the reduced density matrix for the qubits; the 
dephasing terms are extracted from the off-diagonal el- 
ements which are also called the coherence terms of the 
system. The charge fluctuations approacb2ii2^i22 consists 
in investigating the effect of the transit of many electrons 
through the arms. This transport is stochastic and gives 
rise to fluctuating potentials in the arms also seen by 



the DQDs through the coupHng. As a consequence the 
DQDs will exhibit a fluctuating phase which will lead to 
dephasing. 



A. Entanglement approach 

To get the information concerning the dephasing rates, 
we have to write the reduced density matrix of the two 
DQDs at time t, after the transit of N electrons. For this, 
we first consider the transit of one electron (one scatter- 
ing event) in a time t -^ t, t being the traversal time of 
one electron in the MZI. We assume that all consecutive 
scattering events are independent. As a consequence the 
entanglement generated by the measurement will be the 
product of the entanglement generated by each indepen- 
dent electron that transfers through the MZI during the 
measurement time, t = N ■ t. Initially the MZI and the 
two DQDs are completely disentangled: 

= |*dct) ® (alt) + /3|;» ® (7lt> + m). (46) 

where a, /3, 7, 5 are the complex numbers defined in Sec- 
tion V and IvPdct) is the state of the detector before the 
experiment starts. When one electron transits through 
the arms of the MZI, it acquires information about the 
state of the DQDs. As a consequence the state of the 
detector and the state of the qubits become entangled. 
Assuming that the electron can go out in lead 3 and 4 of 
the MZI, the entangled state of the MZI and the qubits 
at time r reads)^ 

|*(r)) = a^e-^^-'\sll\A) + Sllm®\^) 






(47) 



1 3) and |4) are the output detector states and Sij are 
the amplitudes of probabilities to go from input j to out- 
put i. These amplitudes depend on the state of the DQDs 
through the phases xi a-nd X2 (see Eq. ([8]) for the expres- 
sion of 5*31 for instance). The phase factor depending on 
the qubit bias e describes the free evolution of the qubits 
during the time t. The reduced density matrix for the 
DQDs at time t is given by 

p'^\t) = (3|vl/(r))(vl/(r)|3) + (4|v&(r))(vl/(r)|4). (48) 

The elements of this reduced density matrix can be ex- 
pressed in the form)^ 



P?(-) 



P^(0)e- 



•ijT A. 



(49) 



For simplicity the density matrix is written in the com- 
putational basis lltt), Itl), lit), lU)}- pfj{'^) are the el- 
ements of the density matrix p«''(0) = |*9^(0))(*«^(0)|. 
Aij depends on the state of the DQDs through Axi and 



Ax2 and on the reflection and transmission probabilities 
of the left QPC, Rl and T^. The frequency Wy depends 
on e and describes free evolution corresponding to the 
matrix element pij. Then at time t, assuming that the N 
electrons are independent, the reduced density matrix's 
elements are given by: 



ptit)=pfMe- 



*Ai!. 



(50) 



t — Nh/eV where e is the electrical charge and V is the 
applied potential between lead 1 and lead 3. That is to 
say eV/h is the rate at which the electrons are sent into 
the interferometer or in an equivalent manner, it corre- 
sponds to the inverse of the traversal time r introduced 
before. By comparing this expression for the off-diagonal 
elements with the standard expression of the coherence 
terms, pfAt) ~ p? (0)e^"^*J*e~'"''=p''-*^*, we deduce the 
expressions of the dephasing rates: 



dcph.zj 



■loe 



A,: 



N 



(51) 



In the entanglement approach, two or more scattering 
matrices are used, depending on the number of qubits 
in the setup. In our case we deal with 4 scattering ma- 
trices namely depending on the charge state aa' of the 
qubits, 5'^'^ . In contrast, in the charge fluctuations ap- 
proach (see Appendix B), the scattering properties are 
expressed in terms of a reference scattering matrix S 
(scattering matrix in the absence of the qubits) and in 
terms of functional derivatives with regards to potential 
variations. Thus in the entanglement approach, we can 
write: 



*+i*^" 



where SW'' is calculated self-consistentlys 



(52) 



ii. This 

provides the capacitance fraction in the entanglement ap- 
proach. In our case the different scattering matrices differ 
only through phases, i.e Ax = C^/Ci. 

As mentioned in Section IV, we obtain six different 
dephasing rates which reduce to the three given in Eqs. 
(|23II25[) when taking into account the parity assumptions, 
Eqs. dni) and (fTi)) . 



B. Charge fluctuations approach 

This approach investigates the fluctuating potentials in 
the arms of the MZI and in the DQDs, arising from the 
quantum statistical nature of electron transport. In the 
simple case of one two-level system coupled to a detector, 
the dephasing rate is given fa yis, 21, 25^32 



^iHit) 



£/, 



(53) 



where U is the fluctuating potential of the detector and 
(50(i) is the fluctuating phase acquired by the two-level 



system as a consequence of U. Assuming that 6(j){t) is a 
Gaussian phase, the dephasing rate can be written as 



{e'^-PW)^ 



-^{54>itf)u 



(54) 



with {S(j){t)^)u oc Suu{0)t and Suu{0) is the spectral 
density of the voltage fluctuations at zero frequency. 
This spectral density can be expressed in terms of 
a non-equilibrium charge relaxation resistance2ii2^ i?^, 
SuuiO) oc (C^/Ci) Ry. If only a single potential counts 
(spatial variations of the potential can be neglected) , the 
non-equilibrium charge relaxation resistance is given by: 



Ro, — 



h Tr[N,,N*,] ^^J_^^dldU 

e2 (TriV)2 ' 2Tn dU dE' ^ ' 



N \s a. generalized Wigner-Smith matrix (a density-of- 
states matrix) that is determined by the region inside 
the detector where the potential U is non zero. In the 
case of one qubit coupled to a one-channel detector;^ N 
is a 2 X 2 matrix and N12 is its off-diagonal element. 

Adapting this method to the case of two DQDs cou- 
pled to the MZI, we find the dephasing rates given in Eqs. 
(|23ll25p . We assume two fluctuating potentials induced 
by the transfer of the electrons, f/i(i) on the upper arm 
and U2 (t) on the lower arm. Through the coupling capac- 
itances Ci, both DQDs will acquire a fluctuating phase 
on top of their free evolution phase, 54'i{t) for the upper 
DQD and 54)2{t) for the lower DQD, which will depend 
on the state of the DQDs. Indeed there is no interaction 
when the charges in the DQDs are in the dot far from 
the arms. As in the simple case presented above, the cor- 
relators {6(j)i{t)S(t)j{t))ij are proportional to the voltage 
fluctuations Su-u {^) which are themselves proportional 



to the non-equilibrium charge relaxation resistance R^ 
given by: 



(«j) 



j;fa-)_fe Tr[JVgiV,(f] 



a2 TriV(')Tr7V0') ' 



2Tri 



dxi dEi 

(56) 
The phases xi and X2 are those defined by Eqs. ^ and 
([5|) and the Wigner-Smith matrix can be calculated from 
the scattering matrix given by Eq. ([5]). 

We find the following expressions for the correlators: 



{5Mt)5h{t))u 



-*-1? (t 



(57) 



{Hi{t? 



eVt fC^ 

IT I a 



h2{tf)u = RLTL^-^[^\ .(58) 



These correlators show that the events on the upper arm 
and lower arm are perfectly anticorrelated as expected 
from the conservation of the current. Due to this perfect 
anticorrelation, one finds by calculating the six dephasing 
rates that they reduce to three: 



reo = Y^iHiAtf) 



1 



1 eV fCf, 



(59) 



2^ ^(6Mt)^sMt)y)^^2RLn^(^]m 



r -1 

°° ~ 2t 



L(t)+^02(i))')^=O. 



(61) 



It is interesting to notice that the assumption Axi = 
Ax2 = Ax is equivalent to assume a perfect anticorrela- 
tion between S(j)i{t) and 5(t)2{t). These rates are exactly 
those obtained with the entanglement approach. 
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